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ABSTRACT 

Two submm/mm sources in the Barnard lb (Bl-b) core, Bl-bN and Bl-bS, have been observed with 
the Submilhmeter Array (SMA) and the Submihimeter Telescope (SMT). The 1.1 mm continuum map 
obtained with the SMA reveals that the two sources contain spatially compact components, suggesting 
that they harbor protostars. The N2D+ and N2H+ J=3-2 maps were obtained by combining the SMA 
and SMT data. The N2D"'" map clearly shows two peaks at the continuum positions. The N2II+ map 
also peaks at the continuum positions, but is more dominated by the spatially extended component. 
The N2D+/N2H+ ratio was estimated to be 0.2 at the positions of both Bl-bN and Bl-bS. The 
derived N2D+/N2II+ ratio is comparable to those of the prestellar cores in the late evolutionary 
stage and the class protostars in the early evolutionary stage. Although Bl-bN is bright in N2H+ 
and N2D+, this source was barely seen in H^^CO"'". This implies that the depletion of carbon- 
bearing molecules is significant in Bl-bN. The chemical property suggests that Bl-bN is in the earlier 
evolutionary stage as compared to Bl-bS with the H^^CO^' counterpart. The N2H+ and N2D+ lines 
show that the radial velocities of the two sources are different by ^ 0.9 km s~^. However, the velocity 
pattern along the line through Bl-bN and Bl-bS suggests that these two sources were not formed out 
of a single rotating cloud. It is likely that the Bl-b core consists of two velocity components, each of 
which harbors a very young source. 

Subject headings: Stars: formation - Stars: individual: Bl-bN - Stars: individual: Bl-bS 



1. INTRODUCTION 

The object formed at the beginning of the star- 
formation process, the first hy drostatic core, was the- 
oretically predicted by iLarsonl (1969). Since the first 
hydrostatic core is in the transient stage between star- 
less and protostellar cores, it is a key object to under- 
stand how stars are formed in molecular clouds. In 
the spherical symmetric case, first hydrostatic cores are 
extremely short-lived, making the expected n umber of 
them only few hundredth of Class protostars (jOmukail 
120071: iMasunaga fc Inutsukall2Q00f ). On the other hand, 
recent numerical simulations have s hown that the ro- 
tating first cores can su rvive longer (jSaigo et al. 1 120081 : 
iSaigo &: Tomisakal 120061 ). providing more opportunities 
to observe them. Lately, some candid ates for fir st cores 
have been reported (e.g., [ Chen et ah 2010; E noch et al.l 
[20T0t iPineda et all [201lt IChen et al., 2012): they have 
extremely low luminosity of < 0.1 Lq and are not de- 
tectable at wavelengths shorter than ~ 24 /xm. However, 
the nature of these newly discovered sources have not 
been well understood yet. 

The chemical properties are also useful tools to char- 
acterize the evolutionary stage of prestellar and pro- 
tostellar cores. In the cold (~ 10 K) and dense (~ 
10* cm~^) environment, carbon-bearing s pecies such as 
CO and CS condense onto dust grains ('C aselli et al. I 
[200 2: Tafalla et al. 2002). In contrast, nitrogen-bearing 
species remain in the ga s phas e longer ()Bergin fc Langen 
[l997t iBelloche fc Andrdl200l . It is known that N2H+ 
and N2D+ lines are considered to be good tracers of cold 
dense gas because their abundances increase considerably 
if the CO is depleted. Observations of prestellar cores 



show that the spacial distributions of N2H+ and N2D+ 
well agree with that of dust continuum emission. In addi- 
tion, in the cold and dense environment, the abundance 
of deuterium isotopologues is known to be enhanced sig- 
nificantly as compared to the elemental D/H ratio of ^ 
1.5 X 10-5 (|Oliveirall200l[ ). The N2D+/N2H+ ratio de- 
rive d in a sample of pre s tellar cores ranges from 0.05 to 
0.4 (jCrapsi et al. I [20051: [GctIu et al.ll200l . There is a 
tight correlation between N2D^ /N2H "'" ratio and the CO 
depletion factor (jCrapsi et al. |[2005[ ). Once a protostar 
is formed in the core, the protostar warms up its environ- 
ment and starts evaporating CO into the gas phase. As 
a result, the N2D^/N2H"'" ratio d rops as the protostar 
evolves ([Emprechtinger et al.ll2009i ). These results sug- 
gest that the N2D+/N2H+ ratio can be used as an indica- 
tor of evolution in both prestellar and protostellar cores. 
Especially protostellar cores with high N2D+/N2H+ ra- 
tio are particularly important because they are consid- 
ered to be in a stage shortly after the formation of the 
central source. 

In this study, we present the N2H+ J=3-2 and N2D"*' 
J=3-2 observations with the Submilhmeter Array 
(SMAfl and the Submihimeter Telescope (SMT) toward 
two sources embedded in the Barnard 1 (Bl) dark cloud 
(d=2 50±50 pc; ICernis fc Straiivi [2003t IHirota et al.l 
|2008[) . which is part of the Perseus molecular cloud 
complex. These two sources, labeled as Barnard 1-bN 
(Bl-bN) and Barnard 1-bS (Bl-bS) by iHirano et ahl 

^ The Submillimeter Array (SMA) is a joint project between the 
Smithsonian Astrophysical Observatory and the Academia Sinica 
Institute of Astronomy and Astrophysics and is funded by the 
Smithsonian Institute and the Academia Sinica. 
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(|1999f ). have no counterparts in the Spitzer images 
taken by the "Fr om Molecular Cor es to Planet Form- 
ing Disks" (c2d; lEvans et alj 120031 ) survey in all four 
bands of IRAQ, and 2 4 /xm. 70 /xm band s of MIPS 
(jJorgensen et al. I [2001 iRebuU et al] [2007). Both 
sources reveal the similar SEDs, which show extremel y 
low dust temperatures of < 20 K (|Hirano fc Liull2013l) . 
Recently, iPezzuto et alJ (|2012[ ) detected the far-infrared 
emission from Bl-bN and Bl-bS using the Herschel 
Space Observatory. On the basis of the SEDs, they 
proposed that Bl-bN and Bl-bS are first core candidates. 

2. OBSERVATIONS 
2.1. Millimeter-interferometer Observations — SMA 

SMA observations of the N2H+ J=3-2 line were made 
on 2008 September 3. The array was in the subcompact 
configuration with seven antennas. The primary-beam of 
the 6 m diameter antennas at 279 GHz was measured to 
be ~ 45". The phase tracking center of the N2H''" J—3- 
2 observations was at a(2000) = 3°33'21'.'4 and (5(2000) 
= 31°07'35'.'3. The spectral resolution mode was set to 
101.6 kHz (1024 channels per 104 MHz chunk) for the 
N2H+ J=3-2. The corresponding velocity resolution is 
0.11 km s~^. A pair of nearby quasars, 3C84 and 3C111, 
were used for gain calibration. Uranus was used for band- 
pass and absolute flux calibrations. 

Observations of the N2D"'" J=3-2 line were carried out 
on 2007 September 10**^ and 11'^. The array was in the 
subcompact configuration. The primary-beam size of the 
antenna at 230 GHz was measured to be ~ 54". The 
phase tracking center was at the same position as the 
N2H'*' observations. The frequency resolution was set to 
be 101.6 kHz, which corresponds to a velocity resolution 
of 0.13 km s~^. The N2D''' J=3-2 line was observed 
simuhaneously with the CO J=2-l, ^^CO J=2-l, C^^O 
J=2-l, and 230 GHz con tinuum, which will be presented 
in the separate paper bv iHirano fc Liu I ()2013t ). 

The visibility data were calibrated and edited using 
the MIR software package. The calibrated visibility data 
were Fourier-tran sformed and CLEANed with MIRIAD 
(|Sault et al.llT99l . using natural weighting. The synthe- 
sized beam sizes were 4"53 x 3'.'21 with a position angle 
(P.A.) of 44.0° for the N2H+ J=3-2 and 6"08 x 3"17 
with a P.A. of 49.2° for the N2D+ J=3-2. 

The 1.1 mm continuum data were obtained by aver- 
aging the line-free channels of the N2H+ J=3-2 obser- 
vations. The 1.1 mm continuum map had a resolution 
of 4'.'53 X 3"21, and the rms noise level of ~ 0.35 mJy 
beam~^. 

2.2. Single-dish Observations — SMT 

Observations of the N2H+ J=3-2 and N2D+ J==3-2 
lines toward Bl-b region were undertaken with the Ari- 
zona Radio Observatory sub-millimeter Telescope (SMT) 
in 2008 November and December. In both lines, we 
made a 9 x 9 point map with a grid spacing of 10". 
The map covers an area of 1.5' x 1.5' in right ascen- 
sion and declination centered at a(2000) = 3°33'56" and 
^(2000) = 31°09'34"8. We used the 1.3 mm receiver 
system equipped with sideband separating mixers. The 
main-beam efficiency of the telescope was ~ 0.74. The 
half-power beam widths (HPBW) were 27" and 33" in 
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Fig. 1. — The 1.1 mm dust emission observed with the SMA. The 
contour levels are -3 ,3 ,6 ,12 ,24 ,36 ,48 ,60 ,72 ,90 ,108 ,126 ,144 
times the la sensitivity of 3.5 mJy beam~^. The crosses indicate 
the peaks of 1.1 mm continuum emission and the peak of the spitzer 
source. The SMA primary beam size is shown as the large circle. 

the N2H"'" J=3-2 and N2D"'" J=3-2 observations, respec- 
tively. For the N2H+ J=3-2 observations, we used the 
chirp transform spectrometer (CTS), which provided a 
spectral resolution of 0.04 MHz. On the other hand, the 
N2D"'" J=3-2 spectra were obtained using the filter bank 
spectrometers, the spectral resolution of which was 0.25 
MHz. 

2.3. Combined Interferometric and Single-Dish 
Observations 

We combined the SMT data with the SMA data in 
order to fill the short-spacing information that was not 
sampled by the interferometer. Wc followed the proce- 
dure d escri bed by Xcn_ct al (201 1), Takakuwa et al. I 
(pool and iWilner fc Welch I (| 19941 ). The combination 
of the data and subsequent imaging were made using the 
MIRIAD software package. 

The single dish image cube was Fourier transformed 
into the visibility data. Then the visibility data made 
from the single-dish image cube and those observed with 
the SMA were inverse Fourier transformed simultane- 
ously into the image plane with natural weighting. The 
SMA-I-SMT combined image cube was made with a ve- 
locity interval of 0.3 km s""'^. The synthesized beam sizes 
of the combined maps were 4'.'21 x 3'.'83 with a P.A of 
49.81° for the N2H+ J=3-2 and 6"40 x 4"05 with a P.A 
of 45.75° for the N2D+ J=3-2. The rms noise level of 
the combined N2H"'" data cube was 0.17 Jy beam~^ and 
that of the N2D+ data cube was 0.13 Jy beam^^. 

3. RESULTS 

3.1. Dust Continuum Emission 

Figure [1] shows the 1.1 mm continuum image of the Bl- 
b region after the primary beam correction. Two sources, 
Bl-bN and Bl-bS, separated by ~ 20" are clearly seen. 
The peak flux densities, total integrated fluxes and peak 
positions of Bl-bN and Bl-bS were determined by two 
dimensional gaussian fitting. The peak flux density at 
Bl-bN is measured to be 0.28±0.03 Jy beam~^ and that 
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at Bl-bS is 0.45±0.04 Jy beam"^. The beam decon- 
volved size of Bl-bN is 2'.'2 x 2'/0 (550 AU x 500 AU), 
and that of Bl-bS is 2'/0 x 1'.'2 (500 AU x 300 AU). The 
total fluxes of Bl-bN and Bl-bS are 0.37 Jy and 0.52 
Jy, respectively. The positions of these compact sources 
are a(2000) = 3h33'"21.34% 5(2000) = 31°07'26'.'4 and 
a(2000) = 3^33^21.19% 5(2000) = 31°07'43'.'6, respec- 
tively. The primary beam of the SMA also included the 
mid infrared source located at ~12" west of Bl-bS. Al- 
though this source was detected in all four bands of IRAC 
and the 24 and 70 /im bands in MIPS [J0rgensen et al.l 
I2006D . there is no counterpart of this source in our 1.1 
mm map. The 3 cr upper limit of this source is 10 mJy 
beam~^. 

We estimated the total mass of gas and dust in Bl-bN 
and Bl-bS by assuming the continuum emission to be 
optically thin at 1.1 mm, and a gas-to-dust ratio to be 
100. The masses are estimated as follows: 



where F^, is the flux density obtained from contin- 
uum emission, D is the source distance, and is the 
dust mass opacity derived from K,y=Ko{i'/i'o)^ , where 
Ko is calculated to be 0.01 c m^ g~^ at kq =231 GHz 
(jOssenkopf fc Henningi 119941) . We adopted the P an d 
Tdust from the fits of SEDs in iHirano fc Li^ (|20ll ). 
The value of /3 at Bl-bN is 1.8±0.4 and that at Bl- 
bS is 1.3±0.2. The values of Tdust at Bl-bN and Bl- 
bS are 15.6±2.2 K and 18.6±1.6 K, respectively. The 
total masses of Bl-bN and Bl-bS are calculated to be 
0.31±0.20 Mq and 0.34±0.14 M©, respectively. 

The H2 column densities at Bl-bN and Bl-bS were also 
estimated using the same assumptions of Tdust and kq. 
The resuhant iV(H2) at Bl-bN and Bl-bS are 2.67±0.31 
X 10^^ cm~^ and 4.84±0.72 x 10^^ cm"^, respectively. 

3.2. Molecular Line Maps 

Figure [2] shows the integrated intensity maps of the 
N2H+ J=3-2 and N2D+ J=3-2 lines. The panels at the 
top of the figure are the maps obtained with only the 
SMA data, while the ones at the bottom include the con- 
tribution of the SMT data. The lines are integrated over 
the velocity ranges from Vlsr = 5.4 - 8.7 km s~^ for 
the N2H+ J=3-2 and from Vlsr = 5.7 - 7.8 km for 
the N2D+ J=3-2. The velocities of N2H+ J=3-2 and 
N2D+ J=3-2 are given with respect to the reference rest 
frequencies of 279.51169 GHz and 231.32990 GHz, re- 
spectively. The peaks of N2H+ J=3-2 and N2D+ J=3-2 
emission agree well with those of 1.1 mm continuum emis- 
sion in both SMA and SMA+SMT maps. In contrast, 
no significant N2H+ and N2D+ emission was detected 
at the position of the Spitzer source. From the compar- 
isons between the line profiles of the SMA data and that 
of the combined SMA-I-SMT data, the missing fluxes of 
the SMA observations were estimated to be ~ 80 % and 
~ 67 % for the N2H+ J=3-2 and N2D+ J=3-2, respec- 
tively. The larger missing flux in N2H+ implies that the 
N2H+ J=3-2 emission is more spatially extended than 
the N2D+ J=3-2 emission. 

Figure [3] shows the velocity channel maps of the com- 
bined SMA-hSMT N2H+ J=3-2 data. In the low velocity 
range (Vlsr = 6.0 - 6.6 km s~^), the N2H+ line emis- 



sion is detected around both Bl-bN and Bl-bS regions. 
In this velocity range, the N2H+ emission reveals a cen- 
trally peaked distribution at Bl-bS (especially at Vlsr 
^ 6.3 km s^"'^), while it is more spatially extended around 
Bl-bN. In the high velocity range (Vlsr = 6.9 - 7.8 km 
s~^), the emission is centrally peaked at Bl-bN, while it 
is faint and spatially extended at Bl-bS. 

The velocity channel maps of the combined 
SMA+SMT N2D+ data (Fig. H]) show the same 
trend as the N2H+ maps; the compact components 
at the positions of Bl-bN and Bl-bS appear in the 
velocity ranges of 6.9 - 7.8 km s~^ and 6.0 - 6.6 km s~^, 
respectively. As in the case of the N2H+, the spatially 
extended N2D+ emission around Bl-bN is also seen in 
the low velocity range. 

3.3. Kinematic Structure 

In order to examine the velocity structure of Bl-b, we 
made the position- velocity (P^V) diagrams of the N2H+ 
and N2D+ along the North-South cut through Bl-bN and 
Bl-bS (Fig. [5|). The upper panels of this figure are the 
diagrams obtained with the SMA data alone, while the 
ones in the lower panels include the contribution of the 
SMT data. The upper panels clearly show two velocity 
components; one is the 7.2 km s^^ component at Bl- 
bN and the other is the 6.3 km s^^ component at Bl- 
bS. At the position of Bl-bN, the N2D+ emission has 
a secondary component at Vlsr ~ 6.3 km s~^. This 
component is connected to the 6.3 km component at 
Bl-bS in the SMA-I-SMT combined map shown in the 
lower panels. The SMA-t-SMT combined map of N2H+ 
shows the faint 7.2 km s~^ component at the position of 
Bl-bS, however, this component is not seen in the N2D+ 
map. 

4. ANALYSIS 
4.1. N2H^ and N2D'^ Hyperfine Fitting 

Since N2H+ J=3-2 and N2D+ J=3-2 rotational tran- 
sitions contain numerous hyperfine components, line pa- 
rameters such as central velocity (Vc), line width (Aw), 
optical depth (r), excitation temperature (Tex) need to 
be determined by hyperfine fitting. We assume all the 
hyperfine components to have the same excitation tem- 
perature and line width, with the relative velocity of 
each component fixed to the laboratory value. Since 
Bl-bN has two velocity components (Fig. [5]), it was 
difficult to determine all the parameters from fitting. 
Therefore, we fixed the excitation temperatures to be 
7.6 K for N2H"'" and 6.7 K for N2D"'" , which were deter- 
mined bv lEmprechtinger et al.l (|2009[ ). Since these exci- 
tation temperatures were derived from the N2H+ J=l-0 
and N2D+ J=l-0 single-dish observations with HPBW 
of ~ 27" and 32", respectively, the values are consid- 
ered to be the average values of the Bl-b region. We 
ad opted the N2H"'" J= 3-2 hyperfine components listed 
in iCaselh et alTI (120021) and those of the N2D+ J=3- 
2 listed in IGerin et al.l (|2001[ ). We used the combined 
SMA+SMT data to perform the fitting. The N2H+ data 
were convolved in order to match the beam size of the 
N2D+ data (6"40 x 4'.'05 with a P.A of 45.75°). The 
primary beam of the SMA antennas were also corrected. 

FigureElshows the combined SMA+SMT N2H+ J=3-2 
and N2D+ J=3-2 line profiles and the results of hyper- 
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Fig. 2.— Integrated line maps of (a) SMA N2H+ J=3-2, (b) SMA N2D+ J=3-2, (c) combined SMA+SMT N2H+ J=3-2, and (d) 
combined SMA+SMT N2D+ J=3-2. The crosses indicate the peaks of 1.1 mm continuum emission and the peak of the spitzer source. The 
black oval on the lower left side and the large circle in each panel show the synthesized beam size and SMA primary beam size, respectively. 
Contour levels start from 3(7 by steps of 3a for all figures. The values of la are 0.17 Jy beam"^, 0.13 Jy beam~^, 0.17 Jy beam~^ and 
0.13 Jy beam~^ for (a) to (d), respectively. 



fine fittings at the positions of Bl-bN and Bl-bS. The 
N2D"'" line profiles were fitted well using the optically 
thin assumption, while the N2H~'' lines were obviously 
optically thick. The spectra at Bl-bN has two velocity 
components; one is at V ~ 7.2 km s~^, which is spatially 
compact, and the other is the spatially extended com- 
ponent at V ~ 6.2 km s~^. The line profiles at Bl-bS 
are mostly dominated by the Vlsr = 6.3 km s~^ compo- 
nent. The determined parameters are listed in Table [T] 
and Table [2I In these tables, the main components refer 
to the dominate components at Bl-bN and Bl-bS, which 
correspond to the velocity components centered at ^ 7.2 
km and 6.3 km s~^, respectively. Since the line pro- 
file of N2H+ J=3-2 at the position of Bl-bS shows the 
secondary component at ^ 7.2 km s^^, we included this 
component to the fitting. 

It should be noted that the N2H+ line becomes op- 
tically thin if the higher excitation temperature is as- 
sumed; the optical depth of the primary component of 
Bl-bN becomes less than 1.0 for the Tox > 8.8 K, and 
that of Bl-bS becomes < 1 for the Tex > 10.5 K. On the 
other hand, the optical depth of the N2II+ line becomes 
the largest (~ 10 for Bl-bN and ~ 18 for Bl-bS) if we 
adopt the lower limit of the excitation temperature de- 
rived from the fitting ('^ 6.4 K for Bl-bN and - 7.2 K 



for Bl-bS). 

4.2. N2H+ and N2D+ Column Density 

We estimated the column densit ies of N2H"'" an d N2D^ 
using the method described by iCaselli et al. I ()2002D . 
The results including column densities of N2H+ and 
N2D+, and the iV(N2D+)/iV(N2H+) ratios at Bl-bN and 
Bl-bS are summarized in Table [3l The column densities 
of N2H+ at the positions of Bl-bN and Bl-bS were cal- 
culated to be ~ 1.1 X 10^3 cm-2 and 1.8 x lO^^ cm-^, 
respectively. While column densities of N2D+ at Bl-bN 
and Bl-bS are - 2.9 x 10^2 cm-^ and 2.4 x lO^^ cm-^, 
respectively. The corresponding N2D+/N2H+ ratios was 
~ 0.27 at Bl-bN and - 0.13 at Bl-bS. These values 
are consistent with the N2D^/N2 H"'" ratio of ^ 0.18 in 
the en tire Bl-b region obtained bv lEmprechtinger et al.l 
()2009() . The N2H+ column density becomes lower if we 
adopt higher excitation temperatures (a factor of 6 - 9 if 
Tex = 20 K). However, the N2D+/N2H+ ratio is almost 
constant at ^ 0.3 - 0.4, because both N2H+ and N2D+ 
lines are optically thin at Tcx^ 10 K. On the other hand, 
if we adopt the lower limit of the excitation temperature, 
the column density of N2II+ increases by a factor of ~ 3 
- 4. Since the N2D+ line is assumed to be optically thin, 
the N2D+/N2H+ ratio decreases to ~ 0.06. 

The fractional abundances of N2II+ at Bl-bN and 
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Fig. 3. — Velocity channel maps of the combined SMA+SMT N2H+ J=3-2 data in Bl-bN and Bl-bS. The central velocity at each 
channel is shown at the top left corner. The small black oval at lower right in each panel indicates the synthesized beam size. The crosses 
are the peaks of 1.1 mm continuum emission and the peak of the spitzcr source. Contour levels start from 3cr by steps of 3cr for all figures, 
where la is 0.17 Jy beam~^. 



Bl-bS were calculated to be 3.97±0.62xl0-i" and 
3.71±0.63x 10~^°, respectively, using the N2H"'' column 
densities and H2 column densit i es der ived from the con- 
tinuum data. |j0rgensen et al. I (|2004[ ) derived the abun- 
dance of N2H"'" in prestellar cores and Class sources to 
be ~ 10"* - 10~^°. The N2H+ abundance in these two 
sources is comparable t o those of the samples studied by 
iJorgensen et al. I ()2004[ ). 

5. DISCUSSIONS 
5.1. Deuterium Fractionation 

The N2D+/N2H+ ratio of 0.13 - 0.27 derived in 
the previous section is 4 orders of magnitude larger 
than the D/H ratio of 10~^ in the interstellar space 
(jOliveiral l2003| l. It is known that the deuterium 
fractionation in N2H+ is related to the core tem- 
perat ure and CO depletion (e.g., iDi Francesco et aL~l 
I2007D . In the proton-deuteron exchange reaction, 
H+-HHD^H2D+-hH2 (IMillar et all [T989h . H2D+/H;^ 
tends to be increased at lower temperature. In addi- 
tion, at densities of > 10^ cm~^, the depletion of carbon- 
bearing molecules onto dust grains further increases the 
H2D+/H+ ratio ((Roberts fc Millar ■.2000. ). Since H2D+ 



is the main ingredient of deuteriated molecules, the en- 
hancement of H2D+ increases the abundance of deuteri- 
ated molecules. Therefore, the N2D+/N2H+ ratio be- 
comes an indicator of the physical conditions and molec- 
ular depletion of dense and cold gas, and is useful to 
probe the chemical and dynamical evolution of cores, in- 
cl uding starless c ores a nd class protostars. 

iCrapsi et al. I (|2005j ) observed 31 low-mass starless 
cores and derived the N2D+/N2H+ ratio to be 0.03 
- 0.44. They found a strong correlation between the 
N2D+/N2H+ ratio and CO depletion, suggesting that 
the N2D+/N2II+ ratio increases with the density of the 
core. The higher N2D+/N2H+ values were seen in the 
cores of later evolutionary stage, i.e. those with centrally 
concentrated density distribution and high H2 column 
densities. In the case of the Class protostellar cores 
studied bv lEmprechtinger et all ()2009D . the N2D+/N2H+ 
ratio also ranges from 0.03 to 0.3. Their study showed 
that the N2D+/N2II+ ratio in the prestellar cores sig- 
nificantly decreases as a function of dust temperature, 
suggesting that the ratio decreases as the central star 
evolves and heats the surrounding gas. These results 
suggest that the largest N2D+/N2H+ ratio is expected to 
be in prestellar cores that are close to form protostars or 
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Fig. 4.— Velocity channel maps of the combined SMA+SMT N2D+ J=3-2 data in Bl-bN and Bl-bS. The central velocity at each 
channel is shown at the top left corner. The small black oval at lower right in each panel indicates the synthesized beam size. The crosses 
are the peaks of 1.1 mm continuum emission and the peak of the spitzer source. Contour levels start from 3cr by steps of 3cr for all figures, 
where Icr is 0.13 Jy bcam~^. 



TABLE 1 

Results of N2H+ J=3-2 Hyperfine Fitting 







Main Component 


Second Component 


Core 


Ttot 


Vc (km s"-^) Av (km s"-^) 


Ttot Vc (km s--^) Av (km s"^) 


Bl-bN 
Bl-bS 


2.00±0.03 
8.00±0.06 


7.19±0.02 0.83±0.05 
6.27±0.02 0.35±0.03 


0.82±0.02 6.18±0.02 0.50±0.05 
O.lOitO.08 7.25±0.01 0.21±0.02 



^ Though the second component is not clear at Bl-bS, two components were used to fit the spectrum 
of N2H+ J=3-2 at Bl-bS because it seems to be a more spatially-extended component centered at ~ 
6.2 km s-l. 



TABLE 2 

Results of N2D+ J=3-2 Hyperfine Fitting 







Main Component 


Second Component 


Core 


P (K) 


Vc (km s--^) Av (km s"^) 


I (K) Vc (km s'^) Av (km s"-^) 


Bl-bN 
Bl-bS 


1.30±0.06 
1.80±0.03 


7.33±0.02 0.94±0.05 
6.43±0.03 0.58±0.05 


0.87±0.08 6.28±0.03 0.69±0.04 



^ The brightness temperature of each hyperfine component is I X (relative intensity of each com- 
ponent). 



TABLE 3 

N2H+, N2D+ AND H2 COLUMN DENSITIES 



Core 


Ar(N2H+) (cm-2) 


Ar(N2D+) (cm-2) 


iV(N2D+)/JV(N2H+) 


Ar(H2) (cm-2) 


iV(N2H+)/Ar(H2) 


Bl-bN 
Bl-bS 


1.06±0.11xl0i3 
1.80±0. 15x10^3 


2. 85±0. 09x10^2 
2.40±0. 09x10^2 


0.27±0.11 
0.13±0.09 


2.67±0.31xl022 
4.84±0.72xl022 


3.97±0.62x 10-10 
3.71±0.63xl0-i'^ 



7 



u 
01 



c 
o 

O 



u 
01 



c 
o 

o 

Q. 




velocity (km/s 



Fig. 5.— P-V diagrams of (a) SMA N2H+ J=3-2, (b) SMA N2D+ J=3-2, (c) Combined SMA+SMT N2H+ J=3-2 and (d) N2D+ J=3-2 
emission along the N-S cut through Bl-bN and Bl-bS. The horizontal hnes indicate the positions of Bl-bN (orange) and Bl-bS (purple). 
Contour levels are from la in steps of 2(t, where Icr are 0.17 and 0.13 Jy beam~^ in the N2H~'' and N2D+ data, respectively. 




5 6 7 

Vlsr [km s-''] 
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Fig. 6.— Line profiles of (a) N2H+ J=3-2 and (b) N2D+ J=3-2 at the position of Bl-bN, (c) N2H+ J=3-2 and (d) N2D+ J=3-2 at 
that of Bl-bS. The dashed curves (orange) indicate individual velocity components, while the solid curves (brown) show the overall fitting 
results. 
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protostellar cores in the earliest evolutionary stage. The 
N2D+/N2H+ ratios derived at Bl-bN and Bl-bS (- 0.2) 
are comparable to those of the evolved prestellar cores 
such as LI 544 and the very young protostars such as 
L1448C, IRAS 03282 and HH211. It is notable that the 
column densities of N2D+ estimated here are the lower 
limits, because we assumed the N2D"'" line to be optically 
thin. Therefore, the N2D+/N2II"'" ratios derived here are 
also the lower limits. 

The high D/H ratios in the Bl- b region ar e also 
shown in other observa tions (|Lis et al . 2002; Roucf f~t al.l 
120051: rMarcelino^et all 12001 . The high D/H ratios are 
consistent with the low temperatu res derived from the 
SEDs of the contin uum emission (jPezzuto et al.l 120121 : 
iHirano fc Liul [20131 ). Since both Bl-bN and Bl-bS con- 
tain compact continuum sources that are clearly detected 
by the interferometer, they are unlikely to be the prestel- 
lar cores. On the basis of the properties of Bl-bN and Bl- 
bS derived from the continuum SEDs such as low temper- 
atures (Tdust < 20 K), low intrinsic luminosities (Ljnt ~ 
0.03 Lq for Bl-bS and < 0.004 L^t, for Bl-b N), and low 
velocity CO outflows (IHirano fc Liu I [20131 ). these two 
sources are considered to be the candidates for the first 
hydrostatic cores. 

5.2. Kinematics 

As shown in Section 13.31 Bl-bN appears at Vlsr = 

7.2 km s~^ and Bl-bS reveals at Vlsr = 6.3 km s~^. In 
other words, two sources with a projected separation of 
^ 5000 AU have a velocity difference of ^ 0.9 km s^^. If 
this velocity difference attributes to the orbital motion of 
gravitationally-bounded objects, the mass of each source 
is estimated to be 0.55 M©. This mass agrees with the 
mass derived from the continuum data within a factor of 
two. If the two sources originate from a single rotating 
cloud, the P-V diagram along the line connecting two 
sources (i.e. in the N-S direction) is expected to show 
a velocity gradient. However, the velocity pattern seen 
in the P-V maps including the short spacing data (Fig. 
[DJc) and Fig. [5][d)) does not show such a velocity gra- 
dient. Instead, these P-V maps reveal a C-shaped pat- 
tern with two velocity components along the line of sight 
of Bl-bN. The spatially extended emission connects Bl- 
bS and the secondary component of Bl-bN, which has 
the same velocity as Bl-bS. In addition, the N2H+ map 
shows an additional 7.2 km s~^ component along the line 
of sight of Bl-bS. Therefore, it is unlikely that the two 
sources were formed in a single rotating cloud through 
the fragmentation. The P-V maps including the short 
spacing data suggest that there are two velocity compo- 
nents along the line of sight of the Bl-b core, and that 
Bl-bN belongs to the 7.2 km s~i cloud and Bl-bS to the 

6.3 km s~^ cloud. 

5.3. Comparison with H^^ CO^ J— 1-0 Emission 

Fig. [Tjshows the P-V diagram oiY(^^CO+ J=l-0 emis- 
sion along the N-S cut through Bl-bN and Bl-bS ob- 
serve d within the NRO 45 ni telescope ( Hirano fc Liu I 
120131) . Two velocity components along the line of sight 
of Bl-b are also seen in the H^^CO"*" map. It is notable 
that the relative intensities of these two velocity com- 
ponents in the H^^CO+ map are considerably different 
from those of the N2H+ and N2D+ maps. As shown in 
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Fig. 7. — P-V diagram of Hl-^CO+ J=l-0 emission along the 
N-S cut through Bl-bN and Bl-bS observed with the NRO 45 
m telescope. The crosses represent the peaks of the 6.3 km s~^ 
component in Bl-bS and the 6.3 and 7.2 km s~^ components in 
Bl-bN. The contour levels are 3, 5, 7, 9, 11, 13, 15, 17, 19, 21 times 
the Icr sensitivity of 0.07 K. 

the previous section, the dense gas around Bl-bN traced 
by the N2H"'" and N2D"'" lines has a velocity of 7.2 km 
s~^. In addition to this 7.2 km s~^ component, there 
is a secondary component at 6.3 km s^^. On the other 
hand, the H^'^CO^ line peaks at 6.3 km s~^ and shows 
no significant enhancement at Vlsr 7.2 km s~^. At 
Bl-bS, both N2H+ and N2D+ maps reveal the 6.3 km 
s^^ component. Although the N2H+ line shows a weak 
secondary peak at ~ 7 km s^^, there is no counterpart of 
that velocity component in the N2D+ line. In the case of 
H^'^CO"'" map, however, it shows double peaks at 6.3 km 
s~^ and 6.9 km s~^. The comparison shows that N2H+ 
and N2D+ emission is partially anti-correlated with the 
H^'^CO+ emission. 

These results can be well explained by the different 
chemical properties between carbon-bearing molecules 
and nitrogen-bearing molecules. Since H^'^CO+ (HCO+) 
is considered to be formed from ^"^CO (CO) through the 
ion-neutron reaction in the gas phase, the abundance of 
H^'^CO"'" is expected to decrease if a significant fraction 
of CO molecules freeze onto the dust grains. In such a 
cold and dense region with CO depletion, the N2H+ and 
N2D+ abundance tend to be enhanced, since CO is the 
main reactant to destroy the N2H+ and N2D+. There- 
fore, the chemical properties of the dense gas surrounding 
Bl-bN with a velocity of 7.2 km s~^ can be explained by 
the scenario of CO depletion. In the case of Bl-bS, the 
H13C0+ coexists with N2H+ and N2D+. This is prob- 
ably because the gas in Bl-bS is heated by the newly 
formed star. As a consequence, the CO in this region 
has been partially come back to the gas phase. The pres- 
ence of CO outflows aroun d this source also support this 
idea ([Hirano fc Liul[20T3l ). Since Hi3C0+ is made from 
the reaction N2H+ + CO ^ ll^^CO+ +N2, the exis- 
tence of N2H"'" and CO makes H^'^CO+. A comparison 
between Bl-bN and Bl-bS suggests that Bl-bN is chem- 
ically youn ger than Bl-bS . The SEDs in IHirano fc Liu I 
(|2013f) and iPezzuto et all (|2012D also show that Bl-bN 
has lower dust temperature. In addition, Bl-bN was not 
detected at 70 /Ltm by Herschel (Pezzuto et al. 2012), 
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while Bl-bS was detected in the same wavelength. These 
results also support the idea that two sources are in 
slightly different evolutionary stages. 

The 6.3 km component along the line of sight of Bl- 
bN is bright in B.^^CO+ while faint in N2H+ and N2D+, 
suggesting that the CO depletion in this velocity compo- 
nent is not significant. On the other hand, the 6.9 km 
s^^ component at Bl-bS is also bright in the H^^CO"'", 
but unclear in the N2H+ and N2D+. Although the LSR 
velocity of this component (6.9 km s~^) is slightly dif- 
ferent from 7.2 km s"^ in the N2H+, the 6.9 km s"^ 
component seen in the H^'^CO^ is probably related to 
the 7.2 km s~^ component observed in the N2H+. The 
chemical properties suggest that the density of the 6.3 
km s~^ component at Bl-bN is higher than that of the 
7.2 km s~^ component at Bl-bS. 

6. CONCLUSIONS 

We have observed the two sources in the Bl-b core, Bl- 
bN and Bl-bS, in 1.1 mm continuum emission with the 
SMA, and N2H+ J=3-2 and N2D+ J=3-2 line emission 
with the SMA and the SMT. The main conclusions are 
the followings: 

1. Two compact sources, Bl-bN and Bl-bS, with a pro- 
jected separation of ~ 20", are shown in the 1.1 
mm continuum map observed with the SMA. Their 
masses were estimated to be ~ 0.3 Mq. 

2. The N2H+ and N2D+ lines show significant enhance- 
ment at the positions of Bl-bN and Bl-bS. 

3. The N2H+ abundance was estimated to be ^ 4 x 
10-1° at both Bl-bN and Bl-bS. The derived N2H+ 



abundance is comparable to those of prestellar cores 
and Class protostars. 

4. The N2D+/N2H+ ratio was estimated to be 0.2 at 
both Bl-bN and Bl-bS. This value is 4 orders of mag- 
nitude higher than that in the interstellar medium. 
The N2D+/N2H+ ratio is comparable to those in the 
late prestellar phase or in the earliest protostellar 
phase. Since these two sources harbor compact com- 
ponents, they are more likely to be in the very early 
stage of protostellar formation. 

5. The radial velocities of Bl-bN and Bl-bS are 
Vlsr=7.2 km s ~^ and Vlsr=6.3 km s respec- 
tively. The C-shaped velocity pattern seen in the 
P-V diagrams of N2H+ and N2D+ does not support 
the idea that two sources are formed by means of 
fragmentation of a single rotating cloud. It is likely 
that there are two velocity components at Vlsr=6.3 
km s -1 and 7.2 km s " along the line of sight of 
Bl-b, each of which harbors a young protostar. 

6. At the position of Bl-bN, the 7.2 km component, 
which is bright in N2H+ and N2D+ , is not clearly seen 
in the H1'^C0+ J=1-0. This suggests that the carbon 
bearing molecules in the dense gas surrounding Bl- 
bN are heavily depleted onto the dust. 

7. On the other hand, the 6.3 km component at Bl- 
bS is observed in N2H+, N2D+ and Hi3C0+ lines. 
The difference in chemical property suggests that Bl- 
bN is in the earlier evolutionary stage as compared 
to Bl-bS. 

We wish to thank all the SMA staff in Hawaii, Cam- 
bridge, and Taipei for their kind help during these obser- 
vations. N. Hirano is supported by NSC grant 99-2112- 
M-001-009-MY3. 
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